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Abstract 

We investigate the Raman spectra in the geometry where both incident and 
scattered photon polarizations are parallel to the z-direction, for a plane- 
chain bilayer coupled via a single-particle tunneling t±. The Raman vertex is 
derived in the tight-binding limit and in the absence of Coulomb screening, 
the Raman intensity can be separated into intraband (oc tj_) and interband 
(oc tj_) transitions. In the small-ij_ limit, the interband part dominates and 
a pseudogap will appear as it does in the conductivity. Coulomb interactions 
bring in a two-particle coupling and result in the breakdown of intra- and 
interband separation. Nevertheless, when t± is small, the Coulomb screening 
(oc has little effect on the intensity to which the unscreened interband 
transitions contribute most. In general, the total Raman spectra are strongly 

dependent on the magnitude of t±. 

PACS numbers: 74.80.Dm, 78.30.-j, 74.25. Jb, 74.72.-h 
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I. INTRODUCTION 



Recently, c-axis properties of the high-T c superconductors have been of increasing in- 
terest both in normal and superconducting states. Among several theoretical works 
a coupled plane-chain model has been proposed |9,10] to study various c-axis properties of 
YBaCuO. Due to a coupling (of magnitude t±) between plane and chain along the c-axis, 
the properties in the c-direction are shown to be quite different as compared to the ones in 
the a-b plane. More importantly, a finite band gap (splitting) between the two renormalized 
bands (arising from the finite coupling) gives a natural explanation of pseudogaps, which are 
mainly observed in the underdoped normal-state copper oxides |TTHl3| and are seen directly 
as a depression of the conductivity at low frequency at low temperatures. 

In the past few years, Raman scattering as a probe to study high T c materials have been 



very successful |fL4]-|27l and have in particular played a significant role in the debate about 
the symmetry of the gap. Owing to the layer behavior of the high T c compounds, most 
experimental and theoretical studies have focused on the a-b plane properties, taking into 
account different symmetrical geometries. In contrast in this paper, we want specially to 
study the c-axis normal-state Raman scattering, based on the coupled plane-chain bilayer 
model. This not only provides another way to test whether this model is valid in high T c 
cuprates, but also gives a path for extending calculations from a simple one band system to 
the more complicated multiband system. In particular, a better theory should take Coulomb 
interactions into account properly for a multiband system. 

Devereaux et al. p8| have recently given a detailed account of Raman scattering in a 
superconducting bilayer with two identical planes. They use the effective mass approxima- 
tion to derive the Raman vertex and have mainly studied the a-b plane properties. In their 
calculations, they have ignored the interband transition which is not so important in the 
a-b plane, but turns out to be the dominant part in some cases for the c-axis. This case is 
discussed in the present paper. 

A summary of what we have done is as follows. We have derived the Raman vertex 
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in the tight-binding limit which is almost exact in a single bilayer limit. We have clearly 
exhibited in the c-axis Raman spectra, the role of the intra- and interband transitions which 
are separable in the absence of the Coulomb screening. The interband transitions contribute 
most in the small-t^ limit. While a mixing term comes in when the Coulomb screening is 
included, the unscreened interband transition still dominates if t± is small. Appendix [A] 
gives the details of how one can separate the response functions needed in the calculation of 
the Raman spectra into intraband and interband contributions. In Appendix ||, we present 
the analogous results for the c-axis Raman spectra in the case of a plane-plane bilayer model. 
In this instance, a simple analytic end expression can be obtained which helps gain physical 
insight into the c-axis Raman process. The total volume has been set to one (Q = 1) in this 
paper. 



II. FORMALISM 

The non-interacting Hamiltonian for the normal state of a coupled plane-chain bilayer 
cuprate is written as 

H = J2clh(k)C k , (1) 

k 

where we have defined the row and column vectors 

ct = (4>4k) ; c k =^ lk j, (2) 

and c lk , c 2k create an electron in the state k in layers 1 and 2 respectively. The energy 
matrix 

/£i(k) t(k z )\ 
h(k) = . (3) 

\t(k z ) 6(k)/ 

Here £i and £ 2 are the (uncoupled) energy bands for the plane and chain respectively and t 
(assumed to be real) is the perpendicular hopping coupling between plane and chain. These 
quantities are defined explicitly later in Eq. (|2"2|). The spin indices are suppressed throughout 
this paper. 



To study the Raman spectra for a multiband system, we use an approach which is suitable 
in the tight-binding limit. The Hamiltonian (Jl|) can be written in the Wannier representation 



as pj 
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i=l Ri,ri 

~ ~k 51 { ct l(R-l)c2(R2)[^R 1+ zd/2,R 2 + 5 Rl _g d / 2 ,R 2 ] + H.C.} , (4) 
Ri,R2 

where Rj are the lattice vectors in a plane, r$ are the displacements to the nearest neighbors 
of Rj, d is the lattice constant in the z direction, and <7j and t± are the nearest-neighbor 
hopping energies within and between the chain and plane layers. The Wannier representation 
is connected to the k-space representation by 

CiiRi) = ^^e^aQc), (5) 
ViV k 

where N is the total number of lattice sites. In the presence of a magnetic vector potential 
A(r), the tight-binding Wannier states are modified by a phase such that |Kj 



Q(Rj) -> Cj(Ri)exp 



-— R» ■ A(Rj 

he 



(6) 



where e is the charge on the electron, % is Planck's constant over 2n, and c is the speed of 
light. The assumption is made that the vector potential is slowly varying over the length 
scale of the crystal lattice and hence A(q) is strongly peaked about q = 0. 
Substituting (||) and @ into @), to second order in A, Eq. @ becomes 
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tf = tf --^p-A+-^A-p-A (7) 

2mc lmc z 



in the k-space, where the vector 



and the tensor 



m v r j dh(k) 



In these formulas m is the bare electron mass. In electronic Raman scattering, we are 
interested in the transition matrix to order A 2 . Therefore, the term A-*p ■ A in ([/]) is used in 
first-order perturbation theory and the term p • A in (0) is used in second-order perturbation 
theory. As a result, Raman scattering probes an "effective density" 

p = E^7(k)a, (io) 

k 

where the vertex matrix has the element 

(i,k|p.e s |Z,k)(Z,k|p.e'|j,k) 



l L 



€j(k) - €i(k) + huj 1 



+ (i,k\p-e%k)(l,k\p-e s \j,k) 



(11) 



ei (k)-e,(k)-^ 

In Eq. (|H|), the summation is over all the uncoupled bands and the dependence of the 
momentum transfer q on 7k is dropped since we assume q <C hp (the Fermi momentum). 
€i(k) is the energy of the i-band, e 7 (e 5 ) is the polarization of incident (scattered) photon, 
and the frequency u 1 {oj s ) denotes the incident (scattered) photon energy. 

For the problem of a plane-chain bilayer, one has two bands. In the case that there is 
a band gap (A) between these two bands which is very small compared to the energies of 
incident and scattered photons, A <C u ! ,u s , the second term of (|TT|) is negligible. The 
Raman vertex matrix in (|T0|) thus reduces to 

7(k) _ m 2 J2cj d2h{k) e s v . (12) 
^ dk fl dk u 

Eq. fll2|) shows how one can study various geometries for the Raman spectra by choosing the 
appropriate incident and scattered photon polarizations. We note that while the vertices 
exhibited in (|12D are somewhat similar to the vertex given by the famous "effective mass 



approximation" pT].|28|, the physical origin is different. 



The differential cross section for Raman scattering is found to be |32 
d 2 a 1 oj s 

— — = -r 2 — r[l + n B (uj)]lmXpp(q -> 0, iu n -> u + i0 + ), (13) 
auoa\l n lj 1 

where nsiyj) = [exp(Pu) — is the Bose distribution function; r = e 2 /mc 2 is the so-called 
Thompson radius, and the Raman density response function 



XpM^n) = ~ [d™^ (T T p q (r)pt(0)), (14) 

which is given in terms of the effective density operator p given in Eq. (|T0|). Taking into 

account the screening effect due to the Coulomb energy 

1 2 

ff c = ^EE £%(qKqaq, (is) 

Z q ij=l 

the Raman response function (JTJj) in flT3|) will be replaced by p^~P^l 

2 2 

X# = Xpp + X] XmUij{<tiXpiP + Yl Xp Pi Uij(<i)Xpjp m U mi(<i)Xpip + ■■• (16) 

ij=l ijml=l 

Both the intralayer and interlayer Coulomb interactions should be included. Here the real 
density on layer i is 



and the response functions xppvXpip are defined analogously to Xpp by Eq. (|14"D . Eq. fllED 
corresponds to RPA which is given diagrammatically in Fig. [I]. 
One can easily reduce the infinite series in (|T6| ) to 



X P 



X% = Xpp + {Xppi,Xpp.){-L-U X r 1 U\ ' 1P \, (18) 

. Xpip ) 

where 1 denotes the 2x2 unit matrix and we have defined the Coulomb interaction matrix 

(Uix U 12 \ 

U=[ (19) 

\U 2 1 U 22 J 

and the real density response function matrix 

(Xpipi Xp\p2 \ 
(20) 
Xp2Pl Xp2p2 I 

In the continuum limit, the intralayer Coulomb interactions U\\ = U22 = ^p- and the 



interlayer Coulomb interactions U12 = U21 = 9||C (<?|| * s the momentum transfer parallel 



2ne 2 g— g||C 

to Orb plane and c is the bilayer spacing ) which all go to infinity in the strong screening 

qy — > limit. In the tight-binding strong-coupling limit, the couplings Uij are themselves 

large. In either cases, one can approximate (|I^) by 

sc -1 / Xpip \ 

Xpp — Xpp ~ ^Xpp\iXpp-i)X I I • (21) 

V Xp2P I 

The minus sign in (^l|) clearly exhibits the screening effect of the Coulomb interaction. 
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III. RESULTS AND DISCUSSIONS 



Eqs. (|13|) and fl2T|) allow one to study various symmetries of the screened Raman intensi- 
ties for two-band systems. In this paper, however, we will focus on the c (or z) axis Raman 
intensities. We first specify the two uncoupled bands and the coupling: 

h 2 
h 2 

t(k z ) = -2t±cos(k z d/2), (22) 

where A corresponds to the band splitting (responsible for the "pseudogap" ) , /i is the chem- 
ical potential, and t± is the plane-chain coupling strength. While the bands £1 and £2 with 
circular and linear Fermi surface are used to simplify the calculation, they give qualitatively 
similar results to those given by more realistic tight-binding bands ||. Since the k z depen- 
dence only comes in through t(k z ) which is assumed to have no k x and k y dependence, the 
only relevant geometry for the c-axis Raman spectra will be (e 5 ,e 7 ) = (z,z). Following 
(0), we define 

7c(**) = , (23) 

W* J 

where j zz = (m/h 2 )d 2 t(k z )/dk 2 = (t±d 2 m/2h 2 ) cos(k z d/2). Eq. fl23|) is the Raman vertex 
used in (p!Q|) for calculating the "c-axis Raman intensities" throughout this paper. 

In Appendix we have shown how one can separate all the response functions of Eq. fl2T|) 
(needed in the calculation of the c-axis Raman spectra) into an intraband and an interband 
contributions. With these separations, one can easily estimate how each susceptibility % 
depends on the coupling strength t±, for both intraband and interband transitions. Of 
most interest, one sees that the unscreened Raman response function xpp has an intraband 
contribution which is proportional tj_ and an interband contribution which is proportional t\. 
Therefore, the interband transitions dominate for the unscreened c-axis Raman intensities 
when t± is small. 
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In Fig. H], we first show the unscreened Raman intensities for small t± = 2meV at different 
temperatures. Fixed parameters A = 20meV and fi = 500meV are used for all the results 
presented in this paper. To include the impurity scattering, we then introduce the scattering 
rates r,- into the Green's functions in (IA6I) 



/ iu n - & + zrisgn(a; n ) -t \ 

G \k,iu; n )= I . (24) 

V -t iw n ~ 6 + iT 2 sgn(uj n ) J 

For simplicity, we will assume that the impurity scattering rates are equal for both bands 

(Ti = T2 = T) and V is linear in temperature (r = 20meV at T = lOOfC), but independent 

of frequency and momentum. As seen in Fig. ^|, the c-axis Raman spectra can exhibit clear 

signatures of a pseudogap, as long as the temperature is not too high. At high temperatures 

for which T > A, the pseudogap feature is washed out by the impurity broadening. This 

is, in particular, shown by the different low-frequency dependence of the intensities (see the 

inset in Fig. 0). 

In Fig. |3], we present the T = 2QK unscreened Raman intensities for different coupling 
strengths. The inset shows that the unscreened c-axis Raman intensity can, in general, be 
separated into an intraband and an interband contribution. For better clarity of presen- 
tation, all intensities have been divided by t\. This leads to a almost constant interband 
contribution plus a intraband contribution proportional to t\. One finds at larger t± that the 
pseudogap feature arising from the interband transitions are overwhelmed by the increasing 
contribution of intraband transitions which tend to dominate the Raman response. 

We next consider the screened Raman intensity based on the full prescription of Eq. ( pT| ) . 



Due to the Coulomb screening, the presence of the second (mixing) term in (pl|) results in the 
breakdown of the intraband and interband separation for the unscreened Raman intensity. 
For comparison with Figs. |2| and ^, we have presented in Figs. |] and [| the screened c-axis 
Raman intensities. Comparing Fig. ^ with Fig. |2], one sees no major difference between 
the unscreened and screened intensities except for different spectral weights of the higher- 
frequency intensity. This result is valid for the small-t^ case and can be established as 
follows. When t± is small, one can ignore xpip 2 ~ t± i n X m (iH) as compared to Xpipn Xp 2 p2 
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~ t° ± (see Appendix |A|). Consequently, (|21|) can be reduced to a more elegant result 

2 (y- ) 2 

^ = X»-E^- (25) 

That is, the screening effects can be separated for each layer. However, the complete sepa- 
ration of intraband and interband contributions is still not possible because the intra- and 
interband parts (both ~ t 2 ± ) of x PPi are equally important (although x PiPi ~ t° ± is dominated 



by the intraband transition). Nevertheless, the mixing (second) term in ( pq) is proportional 
to tj_ and has only a small effect when t± is small. Therefore, the unscreened interband 
transition given by x PP still dominates. For higher values of t±, as used in Figs. |5] and 
|3], the Coulomb screening has a large effect. In fact, in view of the inset in Fig. |5|, the 
Coulomb screening effect tends to cancel against the unscreened interband contribution, 
largely leaving the unscreened intraband transition. 

IV. CONCLUSIONS 

In this paper, we have studied the c-axis Raman spectra such that both the incident 
and scattered photon polarizations are parallel to the z-direction, for a plane-chain bilayer 
two-band cuprate. This coupled plane-chain model |§ has given a natural explanation of 



pseudogaps [|Tl|-[T3[ seen in the AC c-axis optical conductivity. Here we have shown in detail 
how one can also study pseudogaps by doing c-axis Raman scattering. Our calculations are 
based on a tight-binding model for the Raman vertex. Coulomb screening effects are treated 
properly and include both intralayer and interlayer Coulomb interactions. 

When the Coulomb screening is absent, the Raman intensity can be separated into an 
intraband and an interband contribution. The interband contribution (oc dominates in 
the small-t^ limit. However, the presence of the Coulomb screening mixes up this separation 
and the resulting intensities are strongly dependent on the magnitude of t±. Nevertheless, 
when t± is small, the mixing effect arising from Coulomb screening is small and, as a result, 
the interband unscreened transition still dominates. 
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It's obvious that the approach derived in this paper for calculating Raman spectra in a 
two-band system can be easily extended to the case of a multiband system. A extension 
of our calculations of c-axis Raman intensities to the case of a superconducting plane-chain 
bilayer cuprate will be given elsewhere including the effect of impurities. 



ACKNOWLEDGMENTS 

We thank W.A. Atkinson and D. Branch for useful discussions. This work was supported 
by Natural Sciences and Engineering Research Council (NSERC) of Canada and Canadian 
Institute for Advanced Research (CIAR). 



APPENDIX A: INTRABAND AND INTERBAND TRANSITIONS 

In this appendix, we show how one can separate the response functions needed in the 
calculation of the Raman spectra into intraband and interband transitions. Consider a 
general response function 

XAB^iuin) = - fdre^ (T T 4»(T)flt(0)>, (Al) 
Jo 

where here, the operators A, B can be either the effective density p or the real densities p\ 
or p2- These operators can be written in a unified form such that in the limit q — ► 

A = E Ck7(k)C k = £ £ - ; jr: k r, u . (A2) 

k k i,j=l 



In (|A2|) , the matrix 7 is 7 C given in (p3|) for p and 



1 0\ /0 o 1 

71 = I for pi ; 72 = I for p 2 - (A3) 

00/ Vo 1 



Using (|A2|) , the average (• • •) in ([AT]) becomes 



(T T A q (r)St(0)) 



E EE^(k)7^(k / )(T r 4 k (r)c J , k (r)aL, k ,(0) Q , k ,(0)) 

ijml=l k k' 
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E EE7^(k)7^(k')(T TQ , k (0)4 k ,(r))(T T c Jik (r)4 k ,(0)) 

ijml=l k k' 

E E75(k)7^(k)G°(k,-r)G? m (k,r). 

ijml=l k 



(A4) 



The development from the second to the third line in (A3) makes use of the usual Hartree- 
Fock approximation. G ^ are the elements of a 2 x 2 non-interacting Green's function matrix 
defined by 



G (k,r) = -(T r C k (r)C k (0)). 



For the present coupled plane-chain bilayer model, one has (see (Q)) 

' iu n - £1 -t 

-t i(J n - £ 2 , 

Substituting (|A3) into (|A1[) and making use of the Fourier transform 



G x (k, iu r , 



and the orthogonal relation 



1 rP 



13 Jo 



one obtains 



Xab (q - 0, iu n ) = -EEE 1?i (k)7^(k)G° (k, uj^G^Qc, iav + 

P u n i k ijml=l 

Furthermore, using fl23|) , (|A3|) , and the fact that G\ 2 = G^, one finds 



lUO r 



Xpp = ~3 E E Tr [ G o(k, iw n /)7 c G (k, iuv + ^„)7 C ] 

P k 

= X Pl p = ~a E E Tr [ G o(k, ia; n /)7cG r (k ) zu; n , + iw n )7i] 

P "n> k 

X ftPj = Xpj/H = ~3 E E Tr [ G o(k, ia; n /)7iG (k, zav + iwnhj] 

P <"„/ k 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



The trace operator enables one to calculate xppi e ^c in ( |A10 ) in any rotated frame which is 
convenient. We will choose a frame in which the Green's function matrix is diagonal. This 
leads to results that are readily interpretable physically (see later). 



11 



h{k) -> h(k) = W(k)h(k)U(k) 



We first diagonalize the Hamiltonian h in (§) and the Green's function matrix G 1 in 

am 

'e+ 
e_ 
iuo n — e + 
iu n — e_ 

where the unitary matrix is 

U(k)= , " 1 

and the eigenvalues are 



6 + 6 , 
e ± = — = — ± 



+ t 2 . 



In an analogous way, one can rotate the vertices 



7c -» 7c 



f/ f 7cf/ 



7* 



2 * -^-63 



-2i 



(All) 



(A12) 



(A13) 



(A14) 



7i — > 7i 



1 /Ci-e- 1*1 

1*1 e+-£i 



e i - e_ 



and 



72 ^ 72 = ^ f 7 2 [7 



1 ~\t\ 
-1*1 



e i — e 



Go 





/ 7n 


712 \ 


j , 7c = 








\7l2 


-7n/ 




"22 


-"12 


72 = | 








,-"12 


"11 



(A15) 



(A16) 



For convenience, we redefine 

G n 

G 22 

/ «n CK12 

71 = 

\o;i2 a 22 

where one can easily identify these new variables. Substituting ( A17 ) into ( |A10| ), we obtain 



(A17) 
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Xpp = -n 53 53 Til ( G n G 'n + G 22 G 22 ^j + 7i 2 (G11G22 + ^G'n) , 

^ "V k 

Xppi = Xpip = -3 53 53 7u (aiiGnCn - 0(22^22^22) + 7i2«i2 (G u G 22 + G 22 G' U J . 



Xp>2 Xp2P 



o 53 XI ^11 ( a 22^11^n ~ a ll G 2 2G 22 ) - ll20Cl2(GuG 22 + G 2 2G' n 



Xpipi 



— 53 53 ("li^ 11 ^!! + a 22^22C 22 ) + a\ 2 (GxiG' 22 + G 22 G' U 



Xp 2 P2 — ~a 53 Y{ a 22GuG' u + a>l 1 G 2 2G 22 ^ + a^GnGjj + G 22 G' n ^j, 

P "V k 

Xpip 2 = XP2P1 = ~7> 53 53 a n a 22 (CnG^ + G 22 G' 22 ^j — a\ 2 (g u G 22 + G 22 G' x ^j , (A18) 

" "V k 

where G = G(k, iuv) and G' = G(k, iu) n > + zo> n ). Clearly for each \ in ( |A18| ), the first term 
corresponds to an intraband contribution and the second term corresponds to an interband 
contribution. Crudely speaking, when t± is small (recall ( P2| ) and (E3[)), 7n ~ 712 ~ t]_, 
an ~ 0-22 ~ and ai2 ~ These in turn allow one to estimate the ^-dependence 
of each x, for both intraband and interband transitions. 



APPENDIX B: PLANE-PLANE BILAYER MODEL 

In this Appendix, we present the analogous results for the c-axis Raman intensities in a 
coupled plane-plane bilayer. Due to the identity of these two layers, analytical results are 
available and hence gives more insight into the numerical results of the plane-chain model. 
Similar to Eq. (^), we use the simple bands 

h 2 

Uk x ,k y ) = —(k 2 x + k 2 y )-fi + A 

Uk x ,k y ) = ^-{kl + k 2 y )-^ (Bl) 

for the two planes which are separated by a band gap A (pseudogap). The coupling between 
these two planes is the same as given in (p^). In the case of small t±, one can approximate 
e+ = 6 and e_ = & (see (|AT3|) ). 

With the above simplification, the unscreened c-axis Raman intensity given by the imag- 
inary part of Xpp in ( |A18|) can easily be worked out to be proportional to 
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t 2 t uj 



6t 2 ± 4r / 2T 2T 

A 2 (hcu) 2 + AT 2 + {(hco + A) 2 + 4P + (huj- A) 2 + 4P ; 



(B2) 



where the first term corresponds to the familiar intraband transition (oc t\) and the second 
term corresponds to the interband transition (oc t±). One can see from ( |B2| ) that apart 
from the prefactor t^u), the interband term contains the same factors as the intraband term 
except for a peak shift from uj = to uo = ±A in the Lorentzian form. At zero frequency 
the first Drude like contribution in (|B2|) will go like 6t 2 L /A 2 T while the second (interband) 
term will go instead like 4r/A 2 which is the opposite dependence on T than exhibited by 
the first term. Also for t± —>■ the Drude term is small compared with the second term 
which will then exhibit a pseudogap for hu < A as T — > 0. Using (|B2|) , we plot in Fig. ||the 
temperature-dependent c-axis Raman intensities for the plane-plane bilayer model. All the 
parameters used are the same as those used in Fig. ||. It is seen in Fig. |6] that a pseudogap 
can develop in the c-axis Raman intensities at low frequencies provided the temperature (or 
the impurity damping rate T) is not too high. We remark that, in analogy to the plane-chain 
case, the Coulomb screening has little effect on the small-t^ unscreened intensities given in 
Fig. |. 

The major difference between the results in Fig. for a plane-chain bilayer and those 
in Fig. | for a plane-plane bilayer is that the spectral weight just beyond the pseudogap 
frequency is much higher in the former case than in the latter case. This can be explained 
as follows. For a plane-plane bilayer in which the two bands are separated by a constant 
band gap (A) all over the brillouin zone, the interband transitions will be peaked at u = A. 
In contrast in the plane-chain bilayer, the band gap between the plane and chain bands varies 
at different k points (with minimum value of A) and, as a result, the interband transitions 
are spreaded out for u > A. 
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FIG. 1. Diagrams for the screened Raman response function. Each interaction line comprises 
four parts corresponding to intralayer and interlayer Coulomb interactions. 
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FIG. 2. Unscreened c-axis Raman intensities for t± = 2meV at various temperatures. We 
magnify in the inset the different (^-dependence of the low-frequency intensities. 
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FIG. 3. Unscreened c-axis Raman intensities at T = 20K with different values of plane-chain 
coupling t±. For easier comparison, all intensities are divided by t\. The inset shows the sepa- 
ration into intraband and interband transitions of the unscreened c-axis Raman intensity (for the 
t± = lOmeV case). 
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FIG. 4. Screened c-axis Raman intensities for t± = 2meV at various temperatures (to be 
compared to Fig. pi). 
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FIG. 5. Screened c-axis Raman intensities at T = 20K for various values of t± (to be compared 
to Fig. ^). The inset (for t± = lOmeV) shows separably the three different contributions, namely 
the unscreened intraband, unscreened interband, and screening effect. 
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FIG. 6. Unscreened c-axis Raman intensities in a coupled plane-plane bilayer model for 
t± = 2meV at various temperatures. The inset shows the different low-frequency w-dependence of 
the intensities. This is to be compared with Fig. [2] which applies to the plane-chain bilayer. 
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